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ABSTRACT 

We re-examine the method of deprojection of proper motions, which has been used 
for finding the velocity eUipsoid of stars in the nearby Galaxy. This method is only 
legitimate if the lines of sight to the individual stars are uncorrelated with the stars' 
velocities. Very simple models are used to show that spurious results similar to ones 
recently reported are obtained when velocity dispersion decreases with galactocentric 
radius in the expected way. A scheme that compensates for this bias is proposed. 
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1 INTRODUCTION 



iDehnen fc BinnevI l| 19981 . hereafter DB98) introduced a 
method for deprojecting proper-motion data, which allowed 
them to explore the velo city distrib ution of nearby stars in 
the Hipparcos catalogue (|ESA||1997| 1. without knowing their 
radial velocities. This works by taking a weighted ensem- 
ble average of the proper motions of stars found in different 
parts of the sky, under the assumption that the velocity 
distribution is uncorrelated with position on the sky. This 
assumption was legitimate in the case of the sample studied 
by DB98 because all its stars lay within ~ 100 pc of the Sun, 
so it was reasonable to approximate the full phase space dis- 
tribution function by the velocity-space distribution at the 
Sun: /(x,v) :^/(xfl,v). 

Recently iFuchs et all (|2009l . hereafter F09) used the 
DB98 technique to study a sample of stars taken from th e 
Sloan Digital Sky Survey (SDSS: lAbazaiian et all 120091 ). 
This data set contains stars that extend up to ~ 800 pc 
above the plane and span a range of galactocentric radii 
~ 2kpc wide. Since the velocity dispersion of stars varies 
with both radius and distance from the plane, the validity 
of the assumption that the velocity distribution is uncorre- 
lated with sky position is questionable for this spatially ex- 
tended sample. In this paper we demonstrate that applying 
the DB98 technique leads to erroneous results, particularly 
with regard to the tilt of the velocity ellipsoid with respect 
to the Galactic plane. 

In Section [5] we briefly explain the DB98 method, and 
in Section [3] we demonstrate that for a sample like that of 
F09 it gives a biased estimate of the tilt of the velocity ellip- 
soid. Section 13.11 explains the origin of this bias physically. 
Section [3.21 proposes a technique for removing the bias. In 



SectionUwe discuss biases in the DB98 technique more gen- 
erally. 



2 DEPROJECTION 

The deprojection equations are stated and explained by 
DB98, and written out in full by F09. We repeat them here 
for clarity. 

We work in a Cartesian coordinate system, centred on 
the Sun, in which the a;-axis points towards the Galactic 
centre, the y-axis points in the direction of Galactic rotation, 
and the z-axis points towards the north Galactic pole. Given 
a star moving with heliocentric velocity v = (f/, V, W), the 
observed proper-motion velocity is 

p = V-U||S, (1) 

where s is the unit vector pointing from the Sun to the star, 
and is the component of v parallel to s. This can be 
written in matrix form as 



p = A ■ V, where Aij — 5ij — SiSj 



(2) 



The velocity ellipsoid is defined by both the mean veloc- 
ity and the velocity dispersion. To determine the velocity 
dispersion tensor we use the equation 



PiPi = ^ AikVk Ajivi 

k,l 

= B^jkivkvi, 



where 



Bijki = -^{A-ik Aji + Ail Ajk) 



(3) 



(4) 
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is the part of AA that is symmetric in its last pair of indices. 
We are interested in situations in which we know p and 
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Figure 1. Position of the 8 counting volumes, shaded black or 
grey alternately as we look further from the plane, that make up 
a cone with its vertex at the Sun, and its axis perpendicular to 
the Galactic plane. 



§ (and therefore A and B) but do not know v. It is clear 
from the definition of p (equation [IJ that in this case we 
cannot find v for an individual star because we do not know 
v\\. This is reflected in the fact that A is singular. 

We average equations and @ over a sample of stars. 
If the velocities v of these stars are uncorrelated with their 
sky positions s, they will be uncorrelated with A and B, and 
the expectation value of a product such as A ■ v will equal 
the expectation of A times the expectation of v. That is, 
when the velocities are not correlated with § 



{p) = {A.v) = {A).(v). 



(5) 



Provided the stars are sufficiently widely spread on the sky, 
the matrix (A) is not singular, so we can write 

{v) = {A)-^.{p). (6) 

Similarly, 

(vv) = (B)-^ . (pp). (7) 
3 TESTS 

In this section we demonstrate the danger of using the DB98 
technique when the key assumption of uncorrelated v and § 
does not hold. We do this by considering a simplified form of 
the problem addressed by F09, which was to find the velocity 
ellipsoid of stars in the SDSS survey volume. Fig. [T] shows 
our idealisation of the F09 counting volumes - we take them 
to be slices of a cone with the Sun as its apex. In our usual 
Cartesian coordinate system centred on the Sun, the cone 
is defined by VX^+T^ < 0.8Z and Z < 800 pc. It is split 
into eight counting volumes that are each 100 pc thick (cf 
Fig. 5 of F09). 

F09 validated their use of the DB98 method by drawing 
a velocity for every star in their sample from a Schwarzschild 
distribution with a velocity ellipsoid that was everywhere 
aligned with the X, Y and Z axes, and had constant axis 
lengths au, ov and aw- This velocity distribution does not 
vary with position, so v and s will be uncorrelated. In reality 
the velocity ellipsoid will vary from point to point, both in 
the orientation of its principal axes, and in the lengths of 
these ajces. The lengths of these axes are expected to vary 
with galactocentric radius R roughly as cr oc exp{—R/Ra), 
where Ra is of order twice the disc's scale length Rd (e.g. 



iBinnev fc Tremaind 120081 ). In the Milky Way, Rd ~ 2.5 kpc 
(e.g. Juric et aL[2008h r^ Ra ~ 5 kpc. 

To illustrate the difficulty we adopt the distribution 
function 



/ oc exp 



VR 



aR(R,z) 



+ 



— Vc{R) — {v^{R,z)) 
o-^{R, z) 



+ 



a^(R, z) 



(8) 
(9) 



where again [R, (j), z) are cylindrical coordinates centred 
on the Galactic Centre, Vc{R) is the circular speed, and 
{v^{R,z)) is the asymmetric drift. The velocity ellipsoid for 
this distribution function is aligned with the cylindrical co- 
ordinate axes. We assume that we can correctly compen- 
sate for the circular vel ocity using Oort's constants (e.g. 
iFeast fc Whitelockl Il997l ). In all cases we take constant 
{v^{R,z)) = — 26kms~^. In practice (v^) varies with cr, 
but this makes virtually no difference to these results, so we 
ignore it for simplicity. 

We consider the following three forms for cr: 

(i) Constant cr = {aR,a^,a^) = (45, 32, 24) kms"\ This 
is nearly the same distribution function used by F09, except 
with the velocity ellipsoid aligned with the cylindrical rather 
than Cartesian axes. 



(ii) Radially varying cr 
(t{R) =(t{Ro) exp[(i?o - R)/Ra] 



(10) 



with Ra — 5 kpc, Ro = 8 kpc, and (t{Ro) taking the same 



value as in case (i) 



(iii) A form that varies both radially and vertically so as 
to provide reasonable fits to the dispersions reported by F09: 

rT{R, z) = (34 + 20z, 23 + 20z, 19 + 30z) km s"^ 

X exp[(i?o -i?)/i?a], (11) 

where R^ — 5 kpc and z is expressed in kpc. 

In each counting volume, we place 100,000 stars drawn 
randomly from a uniform probability distribution over the 
entire volume. We assign each star a velocity randomly cho- 
sen from the distribution function. We then "observe" this 
star, and find its proper motion. This allows us to com- 
pare the values of v and vv we determine from deprojection 
(equations |6]& [7]) to the real values. 

Since we consider everything with respect to the Carte- 
sian axes defined in Section [S] this yields values for (U), 
(UU), {UV}, etc. We can use these values (and the fact that 
the centre of each counting volume lies at X = F = 0) to 
find the velocity dispersions parallel to the cylindrical axes. 



and a. 



4>z ■ 



ajj, (T0 and and the mixed moments (Jr^, 

Note that the mixed moments may be either positive or 
negative. In Figs.[21E]&|3]we plot (Jr^, or^ and a^^, which 
we define by 



(12) 



The two vertex deviations, which describe the orientation 
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Figure 2. Components of the velocity dispersion tensor cr (bot- 
tom) and the velocity eUipsoid tilt angle with respect to the plane, 
a (top) as a function of height above the plane. The figure shows 
the true values from the sample in each counting volume (dotted) 
and the values found by deprojection (solid). The true velocity 
ellipsoid has principal axes aligned with the cylindrical coordi- 
nate directions and axis lengths that are independent of position 
(i.e. cr = const). 



of the velocity ellipsoid with respect to the cylindrical axes, 
can be found from these values as 

* = -iarctan^^^; (13) 

Q = - - arctan ^ (14) 

In each of these cases, the values of {(/), (V) and (W) 
determined from equation ^ are consistent with the true 
values at Rq. 

The lower panels of Figs. (2] [3] & |4] show the values of 
these velocity dispersions and the mixed moments as func- 
tions of distance from the plane for the three distribution 
functions described above: true values are shown by dotted 
lines, while solid lines show values recovered by deprojection. 
We see that deprojection yields reasonably accurate values 
of (T_R, (T0, (Jz, (T_B0 and even when cr varies significantly 
through the counting volumes, so the DB98 procedure is not 
strictly valid. 

However, the value of cthz found by deprojection is ma- 
terially incorrect in all cases, being slightly negative when 
cr does not vary with R, and positive otherwise. The upper 
panels of Figs. [5J [3] & |3] show that these incorrect values of 
fiiz yield values of the tilt angle as large as a ~ —20°. A tilt 
of the long axis of the ellipsoid towards the plane implied 
by a ~ —20° is similar to that seen by F09. Thus our ex- 
periments demonstrate that the F09 tilt could be an artifact 
that arises because the velocity dispersion increases inwards. 
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Figure 3. Similar to Fig. [2] except for cr oc CKp{~R/ Ra), with 
Ra = 5 kpc. Again, dotted lines show the true values, and solid 
lines show those found by deprojection. 
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Figure 4. Similar to Figs. [2] and |3] with cr oc exp(— iJ/iJtr), with 
-Rct = 5 kpc, and with cr varying with z (as can be seen in the 
figure). Again, dotted lines show the true values, and solid lines 
show those found by deprojection. The value of crji^ found here 
is similar to that found in Figure |3] with the angle a being larger 
at large z because the velocity ellipsoid is rounder, due to (Jz 
increasing more than ct^ (c.f. equation [M} . 
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3.1 Physical interpretation 



To understand why a radial gradient in cr leads to an appar- 
ent tilt of the velocity ellipsoid towards the plane, consider a 
simplified case in which there are two fields, both at Galac- 
tic coordinate 6 = 90 — One is at i = and the other is at 
I — 180. The velocity measured by the proper motion, v^, is 
then 



{VR cos 9 + Vz sin 6, at I — 0; 
cos 6 — Vz sin 0, at I — 180. 



(15) 



Since < < 90, both sin 6* and cos 6* are positive. There- 
fore, in the field at Z = 0, is large when vr and Vz have 
the same sign, while in the field at / = 180 it is large when 
they take opposite signs. In the absence of a radial gradient, 
the signature of a tilt towards the plane is therefore larger 
values of at I — than at I — 180. Clearly a radial gra- 
dient in (T mimics this signature in the absence of a tilt. 
Hence if one deprojects under the assumption that there is 
no radial gradient, the algorithm will account for the data 
by reporting a tilt towards the plane. 
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3.2 A workaround 

Given that good sky coverage is essential to the success of 
the DB09 method, one simply cannot assume that the ve- 
locity distribution is the same at the locations of all the 
stars in a sample that reaches out to > 1 kpc from the Sun. 
A remedy that can be considered is to adopt a functional 
form for the radial variation of cr and to use this form to 
correct the observed proper-motion velocities to the values 
they would have had if cr had been independent of position. 
For example, for each star we could calculate a "corrected" 
proper-motion velocity 

p' = (p - A ■ Vcorr) exp[(i? - Ro) / K] (16) 

with an estimate for the true value of the parameter 
Rcr that controls the radial variation of cr (eq. llOp and 
Vcorr ~ V0 -|- {v^}e^ au adjustment for the Solar motion 
and asymmetric drift. Thus defined, p' would be expected 
to average to zero over all directions and to be the proper- 
motion velocity if there were no variation in cr with radius. 

We test this correction by applying it to simulated data 
generated as in case (ii) above. We know the true value of 
Vcorr in this case, so we ignore the relatively minor uncer- 
tainties which are caused by not estimating this correctly. 

The dashed lines in Fig. [S] show the tilt angle a found 
from corrected proper-motion velocities for three values of 
R'^: 4 kpc (short-dashed), 5 kpc (long-dashed) and 6 kpc 
(dot-dashed). In all three cases the corrected data give much 
more accurate results than the uncorrected data (full curve) , 
but the most accurate results are obtained with 7?^ = 6 kpc 
rather than the true value, 5 kpc; with R'^ = 5 kpc we find 
a ~ 1° at the largest values of z because the correction 
does not address the problem that the axes of the velocity 
ellipsoid are aligned with the cylindrical rather than Carte- 
sian axes. A closely related bias is seen when cr is constant 
(Fig. [2]) . Using a value of ii^ = 6 kpc for the correction gives 
a ~ because it wnder-compensates for the bias due to the 
variation in cr, which inadvertently compensates for the bias 
due to the alignment of the velocity ellipsoid's axes. 



Figure 5. Tilt angle with respect to the Galactic plane (a) 
as a function of height above the plane z. The true velocity 
ellipsoid is oriented parallel to the cylindrical axes, with <t oc 
exp(— i?/(5 kpc)). The figure shows the true value of a for the 
sample (dotted), the value found from the proper motions with- 
out applying any correction (solid), and values found applying a 
correction (equation I16|l with R'^ = 4 kpc (short dashed) , 5 kpc 
(long dashed) and 6 kpc (dot dashed). The "true" correction, 
R'^ = 5 kpc does not return the true value of a because it does 
not correct for the fact that the velocity ellipsoid is not aligned 
with the Cartesian axes. 

If we considered the value a well established, we could 
use corrected proper-motion velocities to determine Ra from 
the data. 



4 DISCUSSION 

In this paper we have focused on the tilt of the velocity 
ellipsoid towards the plane, and may have left the reader 
with the impression that, for example, the tilt in the plane or 
the non-mixed terms {gr etc.) are correctly recovered by the 
DB98 technique. While this is true to a good approximation 
in the cases shown here, it is not always true. 

For example, consider the situation described in Sec- 
tion 13.21 in which we need to know the value of Rcr so we 
can compensate for the variation in <t across the counting 
volume. In an approach to the determination of R^ we might 
split the data into two sets, for \l\ < 90deg and \l > 90deg, 
and find a separately for each set - this gives us enough 
information to find R^- However, if the data are split in this 
way, they produce a bias in the values of the non-mixed com- 
ponents of a (as well as the mixed components). This bias 
is in opposite directions for the two data sets, so strongly 
affects the derived value of Ra, but cancels out when the 
two sets are considered together (hence the lack of bias in 
the non- mixed components in Figs. [3] and [3}. 

Similar biases must always be considered when using 
deprojection. In the tests described above, the symmetry of 
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the counting volumes cancelled out the bias in most com- 
ponents of <T, effectively restricting it to cthz. The counting 
volumes of real data sets will not enjoy the high degree of 
symmetry characteristic of our model sets, with the result 
that biases in the values returned by the DB98 method will 
not be confined to (t_rz. 



5 CONCLUSIONS 

In this paper we have demonstrated that the statistical de- 
projection of proper motions cannot be applied straightfor- 
wardly to data spanning a significant volume of the Galaxy. 
This is primarily because the dependence of the velocity dis- 
persion (T on position violates the central assumption of the 
method. 

Using a simple model we have demonstrated that ap- 
plying this method can suggest a large tilt of the velocity 
ellipsoid towards the plane, even if the actual tilt is zero. It 
seems very likely that this effect is responsible for the re- 
markably large tilt, a — —20°, reported by F09. Correcting 
for this effect in the manner discussed in Section [3.21 would 
probably bring this result much closer to the smaller tilt an- 
gles obtained usin g radial velocities (e.g. ISiebert et al]|2008l : 
iBond et al.ll2009D . 

We note, however, that all components of cr other than 
crz were nearly unaffected by this bias in our tests. For 
a realistic survey volume such as that used by F09, these 
biases are likely to be larger than in our tests and in some 
circumstances may materially affect the results. 
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